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SUMMARY 

This  is  an  investigation  of  the  equations  of  nnotion  and 
physical  parameters  involved  in  stabilizing  the  initial  flight 
of  a  vertically  launched  rocket  by  means  of  a  booster  rocket 
pin-connected  below  the  main  rocket.    The  system  is  designed 
to  stabilize  the  flight  in  its  early  stage  before  the  aerodynamic 
control  surfaces  become  effective.    Stability  of  the  system  is 
dependent  on  the  pendulum  action  of  the  booster  rocket. 

The  equations  of  tnotion  were  derived  from  Lagrange's 

generalized  momentum  equation.    The  differential  equations 

thus  obtained  were  not  solved  but  were  tested  for  stability  by 

m«ans  of  Routh's  stability  criteria.     The  ratio  of  the  mass  of 

the  main  rocket,  M^,  to  the  mass  of  the  booster  rocket,  M2» 

Ml  Ml 

was  investigated  for  the  two  values  r-r-  =  1 .5  and  rjr  «  7.75. 

^  Ml  '^ 

The  system  involving  a  mass  ratio- —   =7.75  was  found 

M^ 

to  be  unstable  under  all  conditions.    However,  the  system  in- 

Mi 
volving  a  mass  ratio  —--.  «  1 .5  was  determined  to  be  stable  in 

the  range  1.62  <^  -<  4.54  x  10^°,  whcrey  is  defined  as  the  ratio 

of  the  distance   ^    from  the  center  of  gravity  of  the  booster  M2 

to  the  pin  connecting  the  strut  to  the  main  rocket  M, ,  divided 

by  the  radius  of  gyration,  k2,  of  the  booster  M,.     In  this  range, 

for  any  given  value  of    /^  ,  stability  was  uniquely  determined  by 


one  value  of  the  ratio  a    =    ^  ,  v/herc    /,     is  the  length  of  the 
•trut  fronri  the  main  rocket  Mj^  to  the  booster  rocket  M^.     Thus, 
for  a  given  booster,  stability  is  primarily  a  function  of  the  ratio 
^  ,  and  for  any  given     -^  ,         >^  is  uniquely  determined. 

Although  the  system  was  found  to  be  theoretically  stable 

Ml  . 

for  the  mass  ratio  — -    =  1 .5  ,  the  ratio    <*^    turned  out  to  be  of 

such  great  magnitude  as  to  make  the  system  entirely  impracti- 
cal for  this  particular  nnass  ratio. 


TABLE  OF  CONTENTS 

PART                                       TITLE  PAGE 
Acknowledgements 
Summary 

I  Introduction  1 

II  Derivation  of  Equations  of  Motion  5 

III  Numerical  Investigation  of  Stability  42 

IV  Discussion  of  Results  63 

V  Conclusions  65 
Appendix  79 
Table  of  Symbols  85 
References  86 


-1- 

I.     INTRODUCTION 

When  a  rocket  or  guided  missile  is  initially  launched  there 
is  a  short  period  when  its  flight  is  unstable  due  to  the  fact  that  at 
low  velocities  the  aerodynamic  control  surfaces  are  ineffective. 

Various  procedures  have  been  adopted  to  stabilize  this 
initial  part  of  the  flight.    The  most  obvious  and  most  successful 
of  these  has  been  the  use  of  guide  rails,  or  launching  towers  which 
hold  the  missile  on  its  course  until  it  has  attained  a  velocity  at 
which  the  control  surfaces  become  effective.     But  this  method 
requires  a  bulky  launcher,  and  special  care  is  necessary  to  main- 
tain the  rails  in  proi>er  alignment, 

A  modification  of  the  guide  rail  system  is  the  "short- 
length"  launcher  which  is  really  a  guide  rail  whose  length  has 
been  reduced  to  a  minlxnum  by  greatly  Increasing  the  initial  ac- 
celeration and  thus  the  velocity  of  the  missile.    This  increased 
acceleration  presents  a  great  many  problems  in  component  design 
to  resist  the  tremendous  accelerative  forces. 

During  World  War  II  the  Germans  were  fairly  successful 
with  a  unique  approach  to  the  problem  in  their  well-known  "V-2" 
rocket.     (Reference  I).     The  "V-2"  was  launched  fronn  a  near- 
vertical  position  with  no  external  restrictions.     Instead  it  em- 
ployed four  carbon  vanes  mounted  in  the  jet  stream  and  activated 
by  the  auto- pilot  to  maintain  the  missile  on  its  proper  flight  path. 
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The  great  difficulty  here  was  that  the  vanes  steadily  burned  out 
with  a  consequent  diminishing  of  control. 

The  subject  of  this  thesis  is  a  third  means  of  initial  flight 
stabilization  based  on  the  concept  of  the  action  of  a  pendulum  and 
independent  of  both  external  guides  and  Internal  control  mechan- 
isms. 

If  the  missile  booster  is  constructed  as  a  separate  unit 
and  attached  to  tho  main  rocket  by  a  pin  connection  a  certain  dis- 
tance below  its  own  center  of  gravity  it  can  be  seen  fron^  Figure 
1  that  any  motion  of  the  main  rocket  will  impart  a  proportionate 
motion  to  the  booster.    Specifically,  If  the  main  rocket  rotates 
about  Its  center  of  gravity  In  a  clockwise  direction,  the  connec- 
ting strut  will  cause  the  booster  to  rotate  also  about  its  own  center 
of  gravity  in  a  clockwise  direction.    Thus,  if  the  main  rocket 
turns  off  course  in  a  clockwise  direction  the  booster  immediately 
turns  in  a  clockwise  direction  also,  that  is,  the  booster  tends  to 
re -align  its  thrust  with  the  axis  of  the  main  rocket.     If ,  by  a 
proper  choice  of  moment  arms,  the  booster  can  be  made  to  over- 
shoot this  position  of  thrust  alignment  slightly  the  thrust  will  have 
a  horizontal  component  which  acts  on  the  strut  to  the  main  rocket 
and  causes  the  main  rocket  to  turn  about  its  center  of  gravity  in 
a  counter-clockwise  direction.    Thus,  the  booster  responds  to  an 
error  and  Its  response  tends  to  counteract  the  error.    Such  a  sys- 
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tem  should  produce  a  sinusoidal  flight  path  if  the  proper  moment 
arms  and  damping  forces  were  applied. 

In  this  analysis  frictional  and  aerodynamic  forces  were  ne- 
glected to  siTxiplify  the  equations.     The  only  damping  force  consid- 
ered (other  than  the  inertia  forces  of  the  system)  was  the  jet  damp- 
ing force  or  the  resistance  of  a  jet  streatn  to  rotation.    The  other 
external  forces  considered  were  the  weights  of  the  main  rocket 
and  the  booster,  and  the  thrust  of  the  booster.     (The  motor  of  the 
main  rocket  is  considered  not  to  be  operating  in  this  analysis). 
Taking  into  account  only  these  forces  the  writer  derived  the  equa- 
tions of  motion  of  the  system  in  two  dixriensions  by  means  of  La- 
grange's generalized  momentum  equation.     Because  of  the  com- 
plicated nature  of  these  equations  no  solution  of  them  was  attempted. 
Instead  the  conditions  for  stability  were  investigated  by  means  of 
Routh's  stability  criteria. 

A  hypothetical  systetn  consisting  of  a  five-second  booster 
pin-connected  below  a  "V-2"  rocket  was  first  investigated.    For 


this  system  the  ratio  of  the  mass  of  the  main  rocket  M^,  to  the 
computation  the  lateral  motion  of  the  center  of  gravity  of  the  sys- 


mass  of  the  booster  M>,  is =7.75.     To  reduce  the  work  of 

^  Ma 


tern  was  ««t  equal  to  zero  and  only  the  vertical  motion  of  the 
center  of  gravity  and  the  rotation  of  the  main  rocket  and  the 
booster  were  considered. 
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Since  no  stable  solution  was  found,  a  system  with  mass 

Ml 

ratio s  1.5  was  next  investigated,  again  neglecting  the  lateral 

M2 

motion  of  the  center  of  gravity  of  the  system.    In  this  case  a  range 
of  stability  was  determined  but  it  was  of  such  a  sensitive  nature 
as  to  make  impractical  further  investigation  for  this  mass  ratio 
in  the  general  case  with  lateral  motion  taken  into  account. 
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II.  DERIVATION  OF  THE  EQUATIONS  OF  MOTION  ^ 

In  deriving  the  equations  of  motion,  only  four  external 
forces  are  considered,  namely,  the  weight,  Mj^g,  of  the  main 
rocket,  the  weight,  M2g#  of  the  booster,  the  thrust,  F,  of  the 
booster,  and  the  jet  damping  force,  D,  of  the  booster.  These 
forces  are  indicated  in  Figure  I. 


^1.^ 


Figure  I.  External  Forces  on  Main  Rocket  and  Booster. 
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The  jet  damping  force,  D,  is  a  force  which  resists  the  ro- 
tation of  the  booster  about  its  own  center  of  gravity.     If  the  booster 
is  of  mass  lA^,  then  the  rate 
of  flow  through  the  nozzle  may 
be  written  as  "J^    ,  or  M,.    As 
the  booster  rotates  about  its  cen- 
ter of  gravity  the  tangential  ve- 
locity at  the  nozzle  exit  will  be 
/r  -^^   °*"  ^s  %  .     Therefore 
the  damping  force,  D,  will  be 

equal  to  the  rate  of  change  of  momentum,  or  M^  '  iXfh^  )•    This 
expression  has  the  dimensions  r^    'L*  —  ,  or  -— —  . 

Frictional  forces  and  aerodynaniic  forces  are  neglected  in 
this  analysis . 

Both  the  main  rocket,  Mj,  and  the  booster,  M,.  are  as- 
signed three  degrees  of  freedom,  namely,  lateral,  vertical,  and 
rotational  motion.    Note,  also,  that  only  the  two  dimensional  case 
is  being  treated.     The  coordinates  of  masses  Mj  and  M^  are  indi- 
cated in  Figure  2.     The  connecting  strut  between  M.  and  M,  is  as- 
sumed to  be  of  zero  mass  and  of  infinite  stiffhoss  for  purposes  of 
this  analysis,  and  consequently  the  coordinates  of  the  center  of 
gravity  of  the  system  are  indicated  on  the  straight  line  joining 
Ml  and  M2     in  Figure  2. 


Figure  2.  Schematic  Diagram  of  Coordinate  System  and 
Dimensions  « 
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It  is  assumed  that  the  propulsion  system  of  the  main  rocket, 
Mj,  will  not  be  operating  during  the  initial  stages  of  launching  and 
consequently  M.  is  a  constant  for  this  analysis.  The  booster,  M,, 
will  of  course  be  burning  with  a  consequent  change  of  nr^ass  during 
launching,  and  the  rate  of  change  of  mass  M2.  say  M^*  will  be 
considered,  as  will  be  the  change  of  position  relative  to  M.  and  Mp 
of  the  center  of  gravity  of  the  system. 

The  first  step  of  the  analysis  is  to  eliminate  the  extraneous 
coordinates  q^ ,  q^ ,  q^,  and  q^  in  order  to  write  the  equations  of 
motion  in  terms  of  the  rectilinear  coordinates  of  the  center  of  grav- 
ity of  the  system,  q^,  and  (i^t  and  the  angular  coordinates  q3  and  q^ 
of  Mj  and  M^  about  vertical  axes  through  their  respective  centers 
of  gravity.     (Since  only  the  onset  of  instability  is  of  interest,  the 
analysis  will  be  restricted  to  consideration  of  small  variations  of 
the  angles  qj  and  q^.)    Once  this  step  has  been  completed,  La- 
grange's generalized  momentum  equation  will  be  used  to  deter- 
mine the  equations  of  motion. 

ELIMINATION  OF  SXTRAN20US  COORDINATSS 

The  coordinates  q^ ,  q^ ,  q^  ,  and  qg  are  always  redundant 
because  of  the  rigid  links  in  the  system  and  the  known  position  of 
the  center  of  gravity  of  the  two  masses.  Therefore  these  extran- 
eous coordinates  may  be  eliminated  by  considering  the  constraints 
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elaborated  below. 

The  first  relation  between  the  coordinates  of  the  individual 
centers  of  gravity  and  the  center  of  gravity  of  the  system  may  be 
found  by  taking  moments  about  the  center  of  gravity  of  the  system 
indicated  in  Figure  3. 


PIN 


Figure  3.  Relationship  of  Mj  and  M2  to  the  center  of  gravity  of 
the  system. 


Thus  ^^^("fs-'f^    '  ^^9^f''T^^    "^ 


or  y^^r--f)    =    ^'^f-fr^ 


(1) 


Furthermore  the  horizontal  and  vertical  separations  of 
the  centers  of  gravity  of  M,  and  M,  in  Figure  4  are  fixed  by  the 
length  of  the  line,  ^  ,  joining  the  two,  and  by  the  angle    ,^q    -  // 
which    Jl    nr^akes  with  the  vertical. 

M.3 


Figure  4.  Relationship  of  Mj  and  M,  to  the  connecting  strut. 
In  Figure  4  the  vertical  separation  is 

and  the  horizontal  separation  is 


(2) 


=  -^ 


^ 

f^^ 


t-^; 


(3) 


-11- 

A  fourth  equation  follows  from  a  conaideratlon  of  the 
similar  triangles  indicated  in  Figure  3. 

and,  since  the  denozninators  of  both  sides  of  this  equation  are 
equal  from  Equation  (1),  it  may  be  written 

Equations  (1) »  (2),  (3),  and  (4)  provide  four  linear  rela- 
tions between  the  coordinates  q    .  q^   »  q     »  and  q    ,  and  q    ,  q     , 
q^  ,  and  q^  .    Consequently  the  former  may  be  expressed  in  terms 
of  the  latter  and  thereby  eliminated  from' subsequent  equations. 
Solving  these  equations  the  following  expressions  are  obtained: 

From  Equation  (2)  /^       =       -^    ""^ {fj  ~  ^ )  '*'  ^  (5) 

From  Equation  (3)  /r      ~        ^    ""^  ^^3  - ^J  f  f-r 

or  ^        =      ^  -  ^    ..>  r/,  -  S)  (6) 

Substituting  Equation  (6)  into  Equation  (1) 
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p.  />^^/^^   =   ^'f  ^  ^/^  '^  /^  ^  ^'^  ^5  -  -^y 

Substituting  Zquation  (5)  into  Equation  (4) 


Substituting  Equation  (7)  into  Equation  (6) 

Fronn  liquation  (2) 

Substituting  Equation  (8)  into  this  equation 


(8) 


(9) 
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yT  =r.  -  -^^^  -^  -r/,-/;  ,10) 

These  expressions  for  the  extraneous  coordinates  in  Equa- 
tions (7),  (8),  (9)t  and  (10)  involve,  in  addition  to  the  principal  co- 
ordinates, the  variables    ^j^»  the  distance  between  the  centers  of 
gravity  of  the  two  components ,  and    <^  ,  the  angle  between  the  strut 
^    ,  and  the  line    J^  ,  joining  the  centers  of  gravity  of  the  two 
components.    These  two  variables  will  now  also  be  expressed  in 
terms  of  the  principal  coordinates. 

In  Figure  4  by  the  law  of  cosines 


^  -^X^^X'-^AA  -^^Cfv-p)  (11) 


and  by  the  law  of  sines 


».H  J 


-rp-fs) 


whence      ^-/)    =   /$  "    "/''" f^  -"'^^f^V'^J  <^^) 


-14- 
These  relations  are  unfortunately  complicated.     There- 
fore, to  simplify  them  for  later  use,  the  fact  that  the  variations 
in  q^      and  q^     are  small  will  be  taken  into  account.     Assuming 
these  angles  sufficiently  small  such  that  terms  of  order  q* 
may  be  neglected,  the  following  approximations  can  be  made 


J//7  f.    *■ 


'A  */ 


and 


Ct^j  *■.    :*-  / 


'/: 


From  Equation  (11) 
From  Equation  (12) 

and  since  Jl    ^    J^--^z 


Thus  the  appropriate  approximations  are 


>      (13) 
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Substitution  of  the  results  of  Equation  (13)  Into  Equations 
(7),  (8),  (9)  t  and  (10)  yields  the  following  approximate  relations 
between  the  extraneous  and  the  principal  coordinates 

^  ^  ^   .     ^^    /^- 


^lf/^2 


(14) 


To  derive  the  equations  of  motion  by  the  method  of  La- 
grange the  expression  for  the  kinetic  energy,  T,  must  first  be 
found  for  use  in  the  Lagrangian  equation 


(16) 


(17) 


where  ^  is  the  "generalized  force"  and  q  .    is  ^^  >    (Reference  2). 
The  generalized  force   -»2.  is  employed  rather  than  the  potential 
energy  because  this  is  a  non-conservative  system. 

If  ly   denotes  the  moment  of  inertia  of  the  mass  My    about 
its  own  center  of  gravity,  the  kinetic  energy  of  translation  and 
rotation  for  the  system  may  be  written 
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'^-  T^fs'-  i^'f.^-  i-  ^/^-V^  -/  ^/:r/-^V  (19) 


DERIVATION  OF  EQUATION  OF  MOTION 
FOR  COORDINATE  q^ 

In  applying  Equation  (18)  to  the  direction  of  the  coordinate 
q, ,  the  value  of    -4^    is  first  computed  using  the  kinetic  energy, 
T  ,  of  Equation  (19) . 

^  A  ^^  /"  ^^ 

The  extraneous  coordinates  q     and  q     may  be  elizninated  from 
this  expression  by  differentiating  Equations  (14)  and  (16)  respec- 
tively.   From  Equation  (14) 

and  from  Equation  (16) 

ilL.    =    / 
Therefore 


The  first  term  in  the  Lagrangian  expression  is  then 


-r^V//,-^/J  -  ^^f  (20) 

By  inspection  of  Squation  (19).  the  second  term  of  the  La- 
grangian expression  is 

-dz:  =  o  (21) 

The  generalized  force  ^.  is  most  easily  computed  by  con- 
sidering the  work  done  by  the  external  forces  as  the  coordinate  q. 
is  varied,  all  other  coordinates  remaining  fixed. 

In  general,  Worky  =-^.     A  f^      ,     or    ^-   - 


In  particular,  the  force    <^,   becomes,  referring  to  Figure  2, 
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But  within  the  present  approximations 

hence     ^,    ^     F ^^  t  I>  {ZZ) 

Substituting  the  results  of  Equations  (20),  (21),  and  {Z2^ 
into  Equation  (18),  the  equation  of  motion  in  the  q     direction  be- 
comes 

DERIVATION  OF  EQUATION  OF  MOTION 
FOR  COORDINATE  q^ 

The  equation  of  motion  in  the  q     direction  is  obtained  in 
an  identical  manner . 
From  Equation  (19) 

and,  differentiating  Equation  (15) 


Consequently  (^     =:    / 


Similarly,  differentiating  Equation  (17) 


so  that 

^     =   / 


Therefore 


e  first  term  of  the  Lagrangian  expression  becomes 


Furthermore  -^     =     ^  (25) 

The  generalized  force    ^^^    is,  referring  to  Figure  2, 
^    =  Work    _     f-f^r^/^J^r  ^  ^^^^f^  ~  ^^y^jT^J  a  f^ 


and  since 


"^  f-  r- 


and  <r^j  J^   ~    / 


^i     ~       -'(^,^^z)j     ^    ^    "  ^A  ^^^) 


-20- 
Subatituting  the  results  of  Equations  (24),  (25),  and  (26) 
into  Equation  (18),  the  equation  of  motion  in  the  q^     direction  be- 
comes 

DERIVATION  OF  EQUATION  OF  MOTION 
FOR  COORDINATE  q^ 

The  first  term  of  the  Lagrangian  expression  is  obtained  in 
an  identical  manner  to  that  employed  for  coordinates  q^  and  q^  . 
From  Equation  (19) 

and,  differentiating  Equation  (14) 

so  that 


Similarly,  differentiating  Equation  (16) 


so  that 
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Therefore 


and  the  first  term  of  the  Lagrangian  expression  becomes 


^' 


(f  J  -''^.^^'''■^^-'^'^^^^^^(^■r.-^'^)     US) 


From  Equation  (19) 


:rentiating  Equation  (14) 


(29) 


and,  differentiating  Equation  (14) 


so  that 


>      (30) 


Similarly,  differentiating  Equation  (16) 


fr=f-;^J^f.- 


z-.-^rJ 


•  o  that 


>      (31) 


# 
^ 


But  now  the  values  of  ^/^    and    ^^  cannot  be  evaluated  to 
the  first  order  from  the  linearized  expressions  for  q       and  q^ 
in  Equations  (15)  and  (17)  respectively  because  the  q     term  has 
vanished  in  the  process  of  linearizing  these  equations.     Therefore 
the  differentiation  must  be  carried  out  upon  the  exact  expressions, 
Equations  (10)  and  (8)*  and  the  results  then  linearized. 
From  Equation  (10) 

and  the  time  derivative  is 


f.-r^^A-'-^^^T^A'-^*^^^^^-'''     *"* 


Therefore 


^./Vj'     "^     ^^,  (33) 


The  value  of  ^^  follows  from  Equation  (11) 
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and  upon  differentiation 

which,  by  retaining  only  first  order  terms  may  be  written 


Similarly 


>     (34) 


The  various  derivatives  of  <P  must  be  computed  from 
iilquation  (12) 


^     =    />-/   ^  p   '  -^^^'^' 


-rf  a 


^o     /_ 


^x 


-^---^.-/.;-/--^A-/#  ^^2^-^'^ J 


^^     f^^p^^^  A' 


j£> 


and  approximately 


^/^      ^  ^y      A'-        '^      Jv 


(35) 


Now 
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and  substituting  the  results  of  Equation  (35)  and  linearizing 

Similarly,  fronn  the  results  of  Equations  (35)  and  (13) 

Substituting  Equations  (34),  (35),  (36),  and  (37)  into  Equa- 
tion (33)  and  retaining  only  first  order  terms 


(^' 


To  find      ^/>    the  same  procedure  is  followed.      Thus,  from 
Equation  (8) 


where  &  =  f- -  S 


and  the  time  derivative  is 
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Therefore 


Substituting  Equation*  (34),  (35),  (36),  and  (37)  into  Equation  (39) 
and  retaining  only  first  order  terms 


The  second  tcrnn  of  the  Lagrangian  can  now  be  written  by 
substituting  Equations  (30),  (31),  (3a).  (38),  (39)  and  (41)  into  equa- 
tion (30) 
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The  generalized  force     ^^    is 


Work  J 


Therefore,  referring  to  Figure  5, 


^,  =  -  -^^^  -  '^/^  -  '^-^  ^-'^z^'  -fA.j',)^.^^] 
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Substituting  iCquations  (3),  (9)  •  and  (10)  to  eliininate  the  extraneous 
coordinates  q     ,  q     ,  and  q     ,  and  collecting  terms 


f,-(A-4h'J^ 


f,-U^A)^^'l^ 


Figure  5.  Schematic  diagram  of  external  forces  acting  on  system. 


Carrying  out  the  indicated  differentiation  with  respect  to  q^ 


<^''V--'""'r-)C-^^y^jfj^^-  ^'^-iffU 
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Subatituting  Equations  (34),  (36),  and  (37)  into  this  expression, 
and  using  the  approximations      '"'^  f/  '^7*      "^"^^     '^T^'  ^   '   •  -=^ 
becomes 


Finally,  retaining  only  first  order  terms,  the  expression  for  the 
generalized  force     ^^    reduces  to 


^'    '     -^Tkf^^^'-^'^^^J 


(43) 


Therefore,  by  substituting  the  results  of  Squations  (28), 
(42),  and  (43)  into  Equation  (18),  the  equation  of  nnotion  in  the  q^ 
direction  becomes 
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DCRIVATION  CF  ICOUATION  OF  MOTION 

FOR  COORDINATE  q^ 

The  first  term  of  the  L,agrangian  expression  is  obtained  in 
an  identical  manner  to  that  employed  for  coordinates  q    #  q    « 
and  q^  . 
From  Equation  (19) 

and,  differentiating  Equation  (14) 

so  that 

Similarly,  differentiating  Equation  (16) 

so  that 


Therefore , 
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and  the  first  term  of  the  Lagrangian  expression  becomes 


^  ^  =  ^/  '  ^/^  -  j^^^^  -^'  r-<r;-^/;; 


-  T^Ti^-^-r^f'-^'j-^) 


(45) 


From  Equation  (19) 


and,  differentiating  Equation  (14) 


SO  that 


>     (47) 


Sinnilarly,  differentiating  .Equation  (16) 
so  that 


\      (48) 


^ 


Here  again  the  values  of  j^  and  ~P  cannot  be  evalu- 
ated to  the  first  order  from  the  linearised  expressions  for  q  and 
q  in  Equation  (15)  and  (17)  respectively  because  the  q^    ternn  has 

vanished  in  the  process  of  linearizing  these  equations.    Therefore 
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the  differentiation  must  be  carried  out  upon  the  exact  expressions. 
Equations  (10)  and  (8),  and  the  results  then  linearized. 

The  time  derivative  of  q      ,  defined  by  Equation  (10),  has 
already  been  found  from  Equation  (32) 

Therefore 


The  value  of   -^^    follows  from  Equation  (U) 

and  upon  differentiation 

which,  by  retaining  only  first  order  terms,  may  be  written 


(49) 


^Jy         ^  ^^ 


^fy  A 


^P""^ 
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Similarly 


^  =  4    =    ^^^---^  rrv-r^)J Cf-^ -fs)  L    (50) 


^/ 


The  various  derivatives  of  <S>   must  b«  computed  from 
Equation  (12) 

Thus 

and  approximately 


and,  substituting  the  results  of  Equation  (51)  and  linearizing. 


^y:'^    =    _^   ^..^     ^    --^  (52) 

Similarly,  from  the  results  of  Squatlons  (51)  and  (13), 


I 
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Substituting  Equations  (50),  (31),  (52),  and  (53)  into  Equa- 
tion (49)  and  retaining  only  first  order  terms 

^7      (/^,^Af^)  y^      /^      V*^  A?^^^       ^    'A    0'^ 

or 

To  find  -^^    the  same  procedure  is  followed .    The  time 
derivative  of  q^  ,  defined  by  Equation  (8),  has  already  been  found 
from  Equation  (39): 


Therefore 


^A 


/^^/«^»     'Jf^  Af,t/^^  t/fy  /Jf»'/^t 
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Substituting  iCquationa  (50).  (51),  (52),  and  (53)  Into  Equation  (55) 
and  retaining  only  first  order  terms 

or 

The  second  term  of  the  Lagrangian  can  now  be  written  by 
substituting  Iilquations  (32).  (39).  (47),  (54).  and  (56)  into  £qua. 
tion  (46) 
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The  generaliied  force  «=^^   is 


Therefore,  referring  to  Figure  5, 


Substituting  Equations  (8),  (9)*  ^nd  (10)  into  this  expression  to 
eliminate  the  extraneous  coordinates  q     ,  q     ,  and  q 

and  collecting  terms  and  performing  the  differentiation 

Substituting  Equations  (13),  (34),  (36),  and  (37)  into  this  expression, 
and  using  the  approximations       ^/^r.H:  f^      and        <^<'*J^'i:i     c^  becomes 
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Finally,  retaining  only  first  order  terms,  the  expression  for  the 
generalized  force  ,J^     reduces  to 


^^  -  -^/^^T^^-^y 


(58) 


Therefore,  by  substituting  the  results  of  Equations  (45), 
(5  7),  and  (58)  into  Equation  (18),  the  equation  of  motion  in  the  q^ 
direction  becomes 

-i^'^-r^f^-^'^    =    -V:.^^--^/  (59) 

Thus  Equations  (23),  (27),  (44)  and  (59)  are  the  four  La- 
grangian  equations  of  motion  corresponding  to  the  four  principal 
coordinates  of  the  system. 
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ANALYSIS  OF  THE  EXTERNAL  FORCES  AND  RSDUCTION 
OF  THE  EQUATIONS  TO  DIMENSIONLESS  FORM 

In  carrying  out  the  stability  analysis  the  state  of  vertical 
motion  corresponding  to  the  q     direction  will  be  considered  con- 
stant during  the  perturbation.     That  is,  the  values  of  q^     and  q^ 
will  be  fixed,  and  their  nnagnitudes  will  be  determined  by  the 
time -instant  at  which  the  stability  of  the  systenn  is  investigated. 
Therefore  the  equation  of  motion  in  the  q     direction.  Equation 
(27),  will  not  be  considered  further  in  this  analysis. 

The  thrust  force,  F,  of  the  booster  may  be  expressed  in 
terms  of  the  mass  rate  of  flow,  M^  ,  and  the  specific  impulse  of 
the  fuel,  Ijf,  ,  as 

The  negative  sign  is  necessary  because  M^     is  numerically  a 
negative  quantity  as  it  occurs  in  the  Lagrangian  derivation.     Thus 
the  thrust,  F,  will  become  a  positive  quantity  when  evaluated  nu- 
merically in  the  next  section. 

Also,  it  was  found  in  the  first  part  of  the  derivation  that 
the  jet  damping  force,  D,  is  proportional  to  the  quantity    /^^  %.. 
In  order  to  make  this  a  numerically  positive  quantity,  it  will  be 
written  as 
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The  moment  of  inertia,  I^  ,  can  be  written  as  the  product 
of  the  mass,  M^  ,  and  the  square  of  the  radius  of  gyration,  k^.  , 
that  is 

The  time  derivative  of  the  moment  of  inertia,  I^.  ,  is 

^Z-  ~  <//<•  /^.*  ^     /^    c/f^^  J       .    However  ,  the  second  term  is 

small  compared  to  the  first  term  and  will  be  neglected  in  this 
analysis.     Therefore  the  following  approximation  can  be  made 

Making  these  four  substitutions  in  the  equations  of  motion 
in  the  q     *  q    *  and  q     directions  the  following  equations  are  ob- 
tained from  Equations  (23),  (44),  and  (59),  respectively. 
From  liquation  (23) 


From  Equation  (44) 


(60) 
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From  Equation  (59) 


(61) 


^;.-^-^/5  ^(^^<^je^y.')r: 


^M^u-^}^^'^'  ^'^^--  '^^^^^ 


(62) 


To  write  these  last  three  equations  in  dimensionless  form 
Introduce  the  dimensionless  coordinates  O^    ,  C^  ,  and  Q^  ,  and  the 
characteristic  units  of  time  =  I      »   velocity  5   Ij,,g»  and  acceleration 

=  g- 
Thus 


)>  (63) 


a, 

=  -4- 
V5 

4;  =  ^ 

^•-^ 

'', 

=  ;5 

■>,  =  ^/./i 

4    -  ^-rf. 

<?, 

-;5 

^.  --  ^.^ 

4  -^;/ 
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Making  the  above  substitutions  into  Equation  (60)  and  divid- 
ing the  resulting  equation  by  M,  g,  the  following  dinnensionless  equa- 
tion is  obtained 


I  r^,*^J  177^    ^^   ^,    ^'       /   At.rf^    r.;j    y^^ 

^4^J'-'^)]'^^        =     -  (64) 

Similarly,  substituting  the  values  of  liquation  (63)  into 
Equation  (61)  and  multiplying  the  resulting  equation  by  ^^  ^* 
the  second  dimensionless  equation  is  obtained 

^/^L^  A-^  J(b  ^/'/^-^'^)*jL  A  -  ^^  f,^j^\l(0=o   (65) 
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Finally,  substituting  the  values  of  Equation  (63)  into  Equa- 
tion (62)  and  multiplying  the  resulting  equation  by      ^^^'^      the 
third  dimensionless  equation  is  obtained 


/^'//^  /:ff.]j^^^  ^ ^  ^  ^   S 

;;^quatlons  (64),  (65),  and  (66)  are  the  three  dimensionless 
equations  of  motion  of  the  system  and  contain  the  following  eight 
dimensionless  parameters 

A  -^  -X^  A- 
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III.    NUMERICAL  INVESTIGATION  OF  STABILITY 

Rather  than  attempt  to  find  a  solution  for  iquatlonB  (64), 
(65),  and  (66),  it  was  decided  to  introduce  selected  values  of  the 
first  six  dimensionless  parannetara  and,  treating    -^    and   ■^- 
as  variables,  to  use  Routh's  Stability  Criteria  to  determine 
the  range  of  combinations  of  these  two  variables  in  which  the 
aystexn  would  be  stable . 

To  simplify  the  writing  of  the  determinantal  equation 
necessary  for  the  use  of  the  Routh  criteria,  introduce  the  fol- 
lowing notation  where 

A.  coefficient  of  Q^.  in  Equation  j 

B.  coefficient  of  Q^.  in  Equation  j 
C. .        coefficient  of  Qy  In  Equation  j 

and 

j      =     1  corresponds  to  Equation  (64) 

j      =     2  corresponds  to  Equation  (65) 

j      =     3  corresponds  to  Equation  (66) 

With  this  notation  Equations  (64),  (65),  and  (66)  may  now 
be  written  as 
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C„    C,+  B,.    C.  +  B^,   C,  +  Aj,  Q,  +  B^,  Q^  +  A,,  Q^    =  0 

B.i   C,  +  C^^  G,  +  B,,  Q,  +  A,^  Q^  +  C,^  Q^  +  B^,  C^  -f  A^^  C..,  =  0   |>(67) 

^j,   0,-i-  Cj,  Q,  +  B,,  C,  +  A,,  O,  +  C^,  Q^  +  B^,  6^  +  A^,  Q^  =  0_ 

where,  letting    -^  =  ^   and  ^=  <>L    ,  the  coefficients  in  ii^quation 
(67)  are 
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y^x^M        ^ 


ftr-^/^ 


t/^z.    z/<j 


Tf4^^ 


/^ 


<^-     =     j^f^-^lt-^j^ 
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Inaamuch  as  oscillatory  solutions  are  of  main  interest, 
assume  c/-       '=<^       <S 


Then 


and 


For  the  system  of  homogeneous  differential  equations  (67) 
to  have  other  than  a  trivial  solution,  the  determinant  of  coefficients 
of  C      must  vanish.     That  is 


B„A 


B^,  A  +  A,, 
C,i\  +  Bj«  A  +  A,, 
Cj,  A*  +  B„  X  +  A^; 


By,  A  +  \ 
Cy^  \  +  Bv;i  A  +  A^ 
Q,3  A^  +  B^,  A  +  i\ 


(67a) 


Expanding  the  determinant  of  (67a) 

C„  F  A*^+  (C„  H  +  B„  F)>%  (B„  H  +  C„  J  -  B.^S  +  B„\V)  a'' 
i-  (C„  K  +  B„  J  -  B,^T  +  B„X)A"'+  (C„  R  +  B.,  K  -  B,;,U  +  B„Y)  X 
+  (B„  a  -  B.^V  +  B„Z)  A    =  0 
which  can  be  written 

9^^  "^  ^.  >^  -^  Vj,^  "^  '?,>^  ^  "9^  ><    +  P^-     =0  (68) 


P.  =  C.,  F 

p^  »  C„  H  +  B„  F 

p^  =  B„  H  +  C  J  -  B.^S  +  B„W 

P, 


K  +  B..  J  -  B„  T  +  B,  X 
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p^      =  C„  R  +  B„  K  -  B,,  U  +  B,,  Y 
p       =  B„  R  -  B,^  V  +  B„  Z 


=  C,,Q 


C„w 


=  A,,B,,  +  A^  B^^  -  /v^B,,  -  A^,B^, 

=  B„C^  -  B^,  C3Z 

=  ^^/Cv^  +  E^/  By*  -  A^/C,^  -  B.^By-. 

=  A,,B,,  +  A^^E,,  -A^,B,,.  A,^B,, 


A.-Q«,j 


The  system  will  be  stable  If  the  coefficients  of  Equation  (68) 
and  the  following  terms  of  Routh's  stability  criteria  are  all  of  the 
same  sign  (Reference  3) 


1 .       ^^  ^  /S/> 


3.       ^ 


0^^*.  -^/^>)>^ 


/^/^^ 


>    (69) 
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NUMERICAL  EXy^MPLS  FOR  "V-2"  TYPE  ROCKET 

To  find  a  numerical  solution  to  the  problem  the  "V-2"  rocket 
(weight  approximately  28,000  lb«.)  was  chosen  as  the  main  rocket, 
Mj,  and  a  small  five-second  booster  providing  an  initial  accelera- 
tion of  1-1/2  g  was  selected  as  the  booster,  M2. 

The  dimensionless  parameters  of  the  system  were  then  es- 
timated to  be 

J^  ^    7.75-         -^  ^  z.vs^/o'^ 


-^^  =   o    (at  launching) 

The  stability  of  the  system  was  then  to  be  investigated  for  various 
values  of  the  two  variables  </=  -^  and  /6=  -^-  , 

Because  of  the  inherent  complexity  of  the  constants  of  Equa- 
tion (69)  t  a  vast  amount  of  time  would  be  required  to  investigate 
every  likely  combination  of  the  variables         <<      and      /^      .    How- 
ever, the  writer  optimistically  set  out  to  evaluate  these  constants 
for  various  con^binations  of  a  and  ^  .    Inasnnuch  as  no  combination 
was  found  where  even  the  constants  p     to  p      were  of  the  same  sign, 
much  less  where  the  test  functions  of  Equation  (69)  were  of  the  same 
sign,  it  was  decided  to  reduce  the  complexity  of  the  system  in  order 
to  exp>edite  the  computation. 
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Accordlngly,  O^    and  Q,    v;ere  assumed  to  be  zero,  that  is, 
the  system  was  restricted  to  vertical  and  rotational  movement  of 
the  center  of  gravity.     This  assumption  reduced  the  system  of  equa- 
tions to  two,  namely.  Equations  (65)  and  (66),  and  reduced  the  labor 
involved  in  applying  the  Routh  stability  criteria  substantially.    Once 
a  stable  range  of  a  and^  had  been  determined,  these  values  were 
then  to  have  been  used  in  the  general  case  for  closer  investigation. 

STABILITY  OF  SYSTEM  WITH  TWO  DEGREES  OF  FREEDOM 

If  both  Q,  and  6,  are  assumed  to  be  zero.  Equations  (65) 
and  (66)  may  be  re-written  (using  the  abbreviated  notation  intro- 
duced earlier)  as 

and  assuming       <^ -<<?•*    »  as  in  the  more  general  case,  the  determi- 
nant for  this  case  becomes 


c,^  a'+  B„  a  +  A,,  C^, aV  B^,  ;^  +  A  ^, 

The  expanded  determinant  nnay  then  be  written 


=  0 
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+  (C,z  A^^  +  A„C^,  +  B,^  B,^  -  C^^  A,,  -  A^C„  -  B^^  B^,)  A* 

+  (B,^A^,  +  A,^B^,  -  B^^A^,  -  A^^B^,)A  +  (A„A^,  -  A^^A,,)  a  0       (70) 


Furthermore,  since  the  parameter  -^* —  equals  zero  at  launching, 
the  last  term  of  l^qxiatlon  (70)  becomes 

Therefore  j^lquation  (70)  reduces  to  a  third  order  equation  of  the 
form 


/P//      ^  /9/1     T^   /^   /^     y  /^J     =    o 


(71) 


where 

P,        '  ^31^*9  ^  ^^z^trj  '  ^^z  %^   ~  B^^C^a 

P,       =  ^3..  Ks  +  ^s^^^s  +  B«  B^,  -  C,,  A,^  .  A^,C„  -  B^^  B,^ 

P,      =  B«  \.  -^  ^z®^,  -  »..^3.  -  ^x».3 
and  the  Routh  stability  criteria  require  that  p^  *  p^    •  p    •  P>  *  ^<^<^ 
//f  -  ^ /^3J   all  be  of  the  same  sign.    (Reference  3). 

The  coefficients  of  Equation  (71)  may  now  be  evaluated  by 
substituting  the  parameters  previously  selected  for  the  "V-2"  rocket 
with  five  second  booster.    Collecting  terms  these  coefficients  may 
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be  expressed  as  follows 


iZ    =  -^zv  ■ 


jr.  'y^ '  /c^ 


^=   J-^^^^^.c,)  (73) 


where 


/5   =   ^aT*  ^  i^^*!^^  ^*  (74) 

where 

where 
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Eefore  computing  numerical  valuoa  of  the  coefficients  of 
Equation  (71)  it  will  be  profitable  to  draw  certain  conclusions  from 
an  inspection  of  iilquations  (72)  through  (75). 

Q  is  defined  as  the  quantity  -^  ,  and  /^   is  defined  as  the 
quantity    -^^    ,  where   ^    and  k^     are  inherently  positive  quan- 
tities by  the  physical  nature  of  the  problem.     Therefore,  it  can 
be  concluded  that  for  any  real  solution  of  the  physical  problem, 
both  Q  and/^    must  be  of  the  same  sign,  that  is,  either  both  posi- 
tive, or  both  negative. 

In  Equation  (72)  the  sign  of  p^    Is  Independent  of  /^   and 
takes  the  sign  of  a.    Therefore  numerical  values  of  p^     need  not 
be  deterznined  at  this  stage  of  the  investigation  since  the  sign  of 
p^    will  be  evident  by  Inspection. 

In  Equation  (73)  a,    Is  always  negative,  and  c,    has  roots 
at  ^  =  0.456  and  /^  =  -1.93 0 
Since  both  a,    and  c,    are 
negative  for /^   <    1.93  and 
since   g  is  negative  when  /3 
Is  negative  ,  it  follows 
that  P|     Is  always  posi- 
tive in  this  region. 

Negative   a  yields  negative  p^    fronn  Equation  (72).     Therefore  p^ 
and  p^    will  be  of  opposite  sign  for  all  values  of  /^  <    1.93,  and, 
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by  Routh's  stability  criteria,  stability  is  impossible  for  this  case. 
By  analogous  reasoning  stability  is  impossible  in  the  range 
0<^<   0.456. 

In  Equation  (74)  a^  =  0  when  ^^  =  0,  and  a^     is  negative  for 
all  other  values  of  /^    ,    Furthermore,  b*    has  a  root  at  /^  =2.54x10     , 
and  c,,    has  a  root  at 
^  -  -2.68x10'^  .    Since 
a,  ,  b^    ,  and  c^    are  all 
negative  in  the  range 

it  can  be  concluded  that  stability  is  impossible  in  this  range  by 
reasoning  analogous  to  that  in  the  case  of  p^    . 

In  Equation  (75)  a^   =  0  when  /5  =  0,  and  a,  is  positive  for 
all  other  values  of  /3  , 
Furthermore  b,  has  a 
root  at/^  =  -1.34x10 
and  at     /<^  =  0  .    Also 
c  J  has  roots  at  /^  = 
0.456,    ^  =  -1.92, 
and   /^  =  -2.68x10"^. 
No  conclusions  can  be 
drawn  in  this  case. 


^-^ 


To  sunnmarize  the  conclusions  of  the  preceding  paragraphs. 
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the  possibility  of  stability  need  only  be  investigated  in  the  ranges 
-1.93  <^<    0  and  2.54xl0''*'<^  <  0.456. 

In  evaluating  the  coefficients  p,    ,  p     ,  and  p     the  follow- 
ing procedure  was  adopted.     Selected  values  of  ^    within  the 
positive  range  2.54x10    -C/^  <   0,456  were  used  to  find  the  con- 
stant coefficients  of  Equations  (73),  (74),  and  (75).     Then  the 
positive  roots  of  the  resulting  equations  were  computed  as  func- 
tions of  a.     Negative  roots  were  not  investigated  since   a  must 
be  of  the  same  sign  as  /^  .     These  values  are  tabulated  in  Table 
I. 

Similarly,  values  of  the  coefficients  and  corresponding 
negative  roots  of  Squations  (73),  (74),  and  (75)  were  computed 
and  are  tabulated  In  Table  II. 

The  data  of  Tables  I  and  Table  II  are  plotted  qualitatively 
In  Figure  6  and  Figure  7. 

In  Figure  8  and  Figure  9  are  plotted  the  logarithms  of 
maxinnum  and  minimunn  limits  of  a  for  which  p,    ,  p     ,  and  p 
are  each  positive. 

It  Is  apparent  from  Figure  8  that  there  Is  no  single  value 
of  a  for  which  p,    ,  p     .  and  p^    can  all  be  positive  since   a^,, 
for  positive  p^    greatly  exceeds   a^.,     for  positive  p^      for  all 
values  of  2.54x10'  </3  <     0.456. 

Likewise  ,  In  Figure  9  there  Is  no  single  value  of  a  for 

'I 
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which  p^    »  Pt    •  ^^^  P      «^*^  ^11  b«  negative  since      \a„.^\       for  nega- 
tive p^     greatly  exceeds     lc^„|    for  negative  p     for  all  values  of 
-1.93  </^  <    0. 

Therefore,  it  is  clear  that  .Equations  (65)  and  (66)  do  not 
have  a  stable  solution  for  any  combination  of  a  and/^    for  the 
parameters  chosen. 

STABILITY  CF  SYSTEM  V/ITH  REDUCED  MASS  RATIO 

In  view  of  the  negative  results  obtained  in  the  previous 
case  it  was  decided  to  investigate  the  effect  of  changing  the  par- 
ameter -^     .     The  following  set  of  dimensionless  parameters 
was  then  selected 


=  a.7^  "  /o 


/^zT^ 


JJ/^   ~  /r-"*^ 


^^ 


(at  launching) 


Note  that  -^  was  decreased  from  7.75  to  1,5,  j,*^ 
was  increased  from  1.96x10  to  6.72x10  ,  and  . «_  was  in- 
creased  from  2.48x10      to  8.75x10     .     The  remaining  parameters 
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were  unchanged.     Thus,  this  second  set  of  parameters  describe 
a  system  essentially  the  same  as  that  first  investigated  except 
that  the  mass  ratio    -^-    has  been  markedly  reduced. 

The  coefficients  of  Squation  (71)  are  again  evaluated  by 
substituting  these  new  parameters.  Collecting  terms  the  coef- 
ficients may  be  expressed  as  follows 

p^     ^   ^  ('a.  cC^  i-  ^o)  (76) 

where      a^    =  s:oi^'/o''^ 


.^^^%  ^,; 


(77) 


where    a,      =    :?.^9  >  /o~ 

^'       -       ^'      ~      ^ 

P,      =    ^  ^*  V-  >^^  ^    -   <r^  (78) 

where    a^  =     -  z.zz  * /o    /^ 

r  =     —  fz.ZlM/o'%*'^  Z.'^^./O'^  *  V./^'/o^^     /.C/'/b'"^ 

'         '  /^ 

p^     =     cr^c^^  ^  ^^^    ^    C,  (79) 

where    a_,  =    7.  3z  ^ /o'"^*^ 
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A8  previously  stated,    a  and/^  must  be  of  the  same  sign. 

Again,  the  sign  of  p^    is  independent  ofy   and  takes  the  sign  of  a  . 

In  Equation  (77)  a, 


— ^ 


is  always  positive  .  Also 
c,  has  roots  at  /c^  =  0.68 
and  at  /^   =  -2.84. 

In  Equation  (78) 
a^    =  0  when    /^=  <?,  and 
<fz.      is  negative  for  all 
other  values  of  /^   . 
Furthermore  b;t    has  roots 
at /f=  -5.6x10*"    and  at 
^  =  2.55x10'^.    Also  c^ 
has  a  root  at/^=  -1. 12x10"*"  . 
Since  a^   ,  b^  ,  and  C;,    are  all  negative  in  the  range  ^<^  <  2.55x10 
stability  is  not  possible  in  this  range. 

In  Equation  (79)  a^    =  0  when/?  =  0,  and   <^  is  positive  for 
all  other  values  of  ^    .     Furth- 
ermore \>3   has  roots  at/<^  =  0 
and  at /^  =  -5.62x10      .    Also, 
Cj     has  roots  at  /<^  =  0  .625  , 
^   =  -2.83  and^   =  -LUxlO*"  . 


-/^ 
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It  may  be  concluded  that  the  possibility  of  stability  need 
only  be  investigated  outside  the  range  0<^'^  2,55x10'   . 

As  in  the  previous  case  selected  values  of  both  positive 
and  negative/^  were  substituted  into  Equations  (77) »  (78),  and 
(79) »  to  compute  the  constant  coefficients  of  Squation  (71),  and 
the  resulting  equations  were  then  solved  for  maxinnum  and  min- 
imum allowable  values  of  a.  The  results  are  tabulated  in  Table 
III  and  Table  IV,  and  are  plotted  qualitatively  in  Figure  10  and 
Figure  11. 

Figure  12  and  Figure  13  are  again  logarithmic  plots  o£ 
the  maximum  and  minimum  limits  of  a  versus  the  logarithm 
of  /^  . 

In  Figure  12  It  can  be  seen  that  stability  is  possible  in 
the  range  10   <  /^  <   4.5x10      only  if  a^^,    for  positive  p^       is 
greater  than  <x ^;^      for  positive  p     .     (The  curves  appear  to  co- 
incide.)   However,  an  analysis  of  Table  III  reveals  that  for  all 
values  of  10"  <^  <    0.68  both  c^    and  c,   are  negative  whereas 
a^    is  negative  while  a,    is  positive  in  this  range.     Therefore 

Thus  it  follows  that  a^^.^       is  greater  than  a^^^^    in  this  range 
and  stability  is  innpossible . 

Similarly,  in  Figure  13  stability  is  only  possible  if 
la^^^j      for  negative  p^    is  greater  than    |a^,^l       for  negative  p^  . 
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Howevsr,  an  analysis  of  Table  IV  reveals  that  for  all  values  of 
-2,8-<'/^<r  -10'    both  c^    and  c,    are  positive  while  a^    is  negative 
and  a.j    is  positive  in  this  range  .     Therefore  ,  |*<i^,j=  jj^j  p  ^  J/  ^  \^^*\ 
while    l^^^^l  =  jA|  |> >  |// -  f^^^l         and  thus     |^^,.|      is  greater 
than     |>f'j^,^|    . 

Therefore  it  is  only  necessary  to  investigate  more  close- 
ly the  range  0.68  </^<  4.5x10^   and  -10^ </<^  <   -2.8. 

Returning  to  Equation  (71)  and  multiplying  the  analytic 
expressions  for  the  coefficients  of  p     and  p     the  following  equa- 
tions are  derived  if  only  first  order  terms  are  retained. 


and 
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Factoring     ^     "^^        ^^^^  i:quatlon  (81) 


^     /^x^^/.  ■    /^x)  ^     .^5-  7    /  (82) 

Comparing  Equation  (80)  with  Equation  (82)  it  can  be  seen 
that  both  a^    and  b^    are  multiples  of  a^.   and  b,,   respectively  by  the 
factor    -^     ^a  ^»y    .     For  positive/^    ,  stability  exists  only  if  o^^^^j, 
is  greater  than  q,^^^      .  that  is.  if       ^^     <      ^    ^^^^  ^^^ 
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But  the  factor       -^     ^*-  ^*^      is  nutner- 


ically  negative  so  that  the  condition  for  stability  can  be  written 


Substituting  the  values  of  c^   and  c^    fronci  Equation  (80)  and  Equa- 
tion (82)  this  condition  becomes 

which,  upon  collecting  terms ,  becomes 


Inserting  the  numerical  values  of  the  parameters,  the  second  term 
of  the* denominator  is  of  the  order  10      since   1^'     ~    ^-^^  '  ^'^ 
Therefore  this  term  can  be  neglected  in  connparison  with  the  first 
term  of  the  denominator  and  the  expression  simplifies  to 

^,      '-0^£)f'-4)  (83) 

Solving  this  inequality  with  the  p>arameters  chosen,  the  result  is 
/^    7^       1.62 
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As  an  example,  choose^  s  2.    Inserting  values  from 
Table  III  into  Equation  (78)  and  Equation  (79),  for  ^  =  2  the 
following  values  of  a  are  obtained 


^  -   -^P-^J 


i      jtC  X /a    f  /  -  Z.^^-  /o 


-J 


Thus  p^  changes  from  positive  to  negative  when 
Oa  =s  5.6x10   (1-2.54x10   ),  while  p^  changes  from  negative  to 
positive  when  q^  s  3  ,6x10*^(1-9. SSxlCJ"*),  and  both  p    and  p^  are 
positive  in  the  range  (5  .6xl0*'-5  .52)<  a<  (5  .6xlo''-l.42) .    Re- 
sults with  a  sinnilar  sensitivity  are  obtained  for  larger  values 
of  ^  but  smaller  values  of  a.    (It  can  be  seen  frotn  Figure  12 
that  the  product  of  a  and  ^  is  apprjxinjately  constant  and  equal 
to  10^.) 

There  is  yet  to  be  applied  the  final  Routh  criterion, 
namely  that     /f  -  "^^^  be  positive  also.    Choosing  p     at  its 
smallest  allowable  value,  say  0    ,  and  calling  the  correspond- 
ing value  of  Pj  some  small  positive  value,  say  s,  /yf  -   ^  y  =  ^ . 
rhus  all  the  stability  criteria  are  satisfied  and  the  system  is 
stable.     Therefore,  it  can  be  concluded  that  the  system  is  stable 


-62- 
for  /9  >  1.62  but  that,  for  any  given  value  cA.^  In  this  stable 
range,   a  is  extremely  restricted. 

Finally,  there  remains  the  negative  range  of^  to  be 
investigated,  namely  -10   </^<-2.8.    But,  for  negative  values 
ci  /^  t  ^    must  be  greater  than  -1.62  by  Equation  (83)  In  order 
that  p^  and  p    both  be  negative »    Therefore  stability  Is  not 
possible  for  negative  values  of  /<^ . 
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IV.  DISCUSSION  CF  RESULTS 

Only  two  possible  sets  of  parameters  for  a  launching  sys- 
tem of  this  type  have  been  investigated.     In  the  first  case  where 
the  mass  ratio  was  7.75  it  was  found  that  stability  was  impossible, 
whereas.  In  the  second  case  where  the  mass  ratio  was  1.5  It  was 
found  that  stability  was  theoretically  possible  for  values  of /^  be- 
tween 1.62  and  4.5x10      If  the  effect  of  lateral  motion  of  the  center 
of  gravity  of  the  system  was  neglected.    For  any  given  value  of  ^ 
in  this  range  stability  was  possible  for  only  one  corresponding 
value  of  a . 

Returning  to  the  physical  aspects  of  the  problen^  consider 
now  the  meaning  of  the  variables  a  and  /^  .  /^    Is  defined  as  -^ 
where  ^  Is  the  distance  from  the  center  of  gravity  of  the  booster « 
Mx  ,  to  the  pin  connecting  it  to  the  strut  to  the  main  rocket,  M,  , 
and  k^    is  the  radius  of  gyration  of  the  booster  ,  M^  .     For  any  sys- 
tern  of  this  type  k^   will  always  be  of  the  order  of  one  foot  or 
greater.     Therefore,  ^  = /0 ^z    will  always  be  equal  to  or  greater 
than  /^  .    That  Is,  as  a  first  approximation,    ^   Is  approximately 
numerically  equal  to  /^  .     Thus,  from  a  practical  consideration 
x<^  ,  and  therefore  t^  ,  must   certainly  not  exceed  some  small  num- 
ber of  the  order  of  1  to  10 ,  depending  on  the  size  of  the  booster, 
for  otherwise  the  pin  will  be  located  too  near  the  booster  exhaust. 
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Now  o  is  defined  as^  where  X  is  the  length  of  the  strut 
from  the  pin  In  the  booster  M^    to  the  center  of  gravity  of  the 
main  rocket »  M,  .    Here  again  there  is  a  practical  upper  limit 
for  ^  ,  first,  because  the  strut  was  assumed  to  be  of  zero  mass 
in  the  derivation  of  the  equations  of  motion,  and  second,  because 
the  strut  was  assumed  to  be  of  infinite  rigidity.    A  long  strut 
would  present  structural  problems  in  rigidity  since  it  is  essen- 
tially a  column  acting  under  both  axial  and  bending  loads.    Re- 
stricting ^  to  values  of  the  order  of  a  few  feet  it  is  clear  that  a 
should  not  exceed  100  or  200  at  most,  since  these  values  corres- 
pond to  an  X   of  25  to  100  feet. 

In  view  of  the  foregoing  remarks  consider  the  product  of 
the  two  variables   a  and /^  .     Thus     «^  =  ^^     ^4-      .    Ifk, 
is  of  the  order  of  1  to  10  feet,  say,  and    yfi    is  restricted  to  25  to 
100  feet,  then  it  can  be  seen  that  the  product  ay^   should  not  ex- 
ceed 25  ,    Returning  to  Figure  12  note  that,  in  the  range  1.62  <^ < 
4.54x10     ,   a v^  is  approxixnately  constant  and  equal  to  about  10     . 
Thus,  even  though  theoretically  stable,  this  system  is  of  no  prac- 
tical significance  for  the  mass  ratio  investigated. 
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V.  CONCLUSIONS 

1 «    A  compound  pendulum  system  of  launching  stabilization  is  not 
stable  for  a  mass  ratio  7.75. 

2.  Such  a  system  is  theoretically  stable  for  a  mass  ratio  1.5  if 
the  effect  of  lateral  motion  is  neglected.    However,  even  In  this 
restricted  case  the  system  is  of  no  practical  significance  due  to 
the  excessive  length  requirennents  for  the  strut  from  the  main 
rocket  to  the  booster, 

3.  Since  the  system  is  quite  sensitive  to  the  particular  param- 
eters chosen,  several  more  investigations  of  specific  cases  are 
required  before  any  general  conclusions  can  be  drawn  as  to  the 
practicality  of  this  method  of  stabilization. 
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Figure  6.   Variation  of  p_    »  P^    »  and  p      with  a  for   several  values 
of  positive  ^  .     (  g'   =  7 .  75^ 
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Figure   7.   Variation  of  p     ,  p     ,  and  p      with  a  for  several  values 
of  negative  ^  .     (  ^  =  7<.75) 
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Figure  8.     Minimum  and  maximum  limits  of  a  which  yield  posi- 
tive p      for  positive^    „     (  i!L  =  7.75) 
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Figure  9.  Minimum  and  maximum  limits  of  a  which  yield  nega- 
tive p.     for  negative  /^ .     (    ^~  ^"^S) 
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Figure  10.   Variation  of  jd     .   P^    »  and  p^     with   a  for  several  values 
of  positive  ^  .     {^  -  l-^) 


Figure   11.     Variation  of  p     ,  p     ,  and  p^     with   a  for  several 
values  of  negative    ^  ^     (  ^~  ^•^) 
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Figure  13  o   Minimum  and  maximum  limits  of  a  which  yield  neg 
tive  p      for  negative  ^  „     (  il  =  \^S)f 
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Figure  14„     Schematic  diagram  of  system  showing  inertia 

forces  for  use  in  Newtonian  derivation  of  equa- 
tions of  motion. 
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APPfiNDIX 

NEWTONIAN  DICRIVATION  CF  EQUATIONS  OF  MOTION 

The  equations  of  motion  for  the  two  dimensional  system 
may  be  derived  by  applying  Newton's  Second  Law  of  Motion  pro- 
vided that  proper  consideration  is  given  to  the  inertia  forces  of 
the  system.    In  Figure  14  these  forces  are  indicated  by  broken 
lines.     In  this  derivation  the  mftin  rocket,  M;   ,  and  the  booster » 
M^  ,  will  be  treated  separately  as  free  bodies  and  the  reaction 
of  the  pin  connection  will  be  replaced  in  Figure  14  by  its  two 
components,  f,    in  the  vertical  direction,  and  f^^   in  the  horizon- 
tal direction.     The  pin  is  assumed  to  be  frictlonless  and  there- 
fore no  moment  is  transmitted. 

Treating  the  main  rocket,  M  ,  ,  as  a  free  body  first,  the 
sum  of  the  vertical  forces  acting  on  it  are 

and  therefore 

/  =  /T^  ^  /^^'  .  (1) 

Similarly,  the  sum  of  the  horizontal  forces  acting  on  M,  are 
and  therefore 
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Finally,  the  sum  of  the  moments  acting  on  M,  are 

En=  ^^    ^  ^^--^  -  /"-^--^ 

or 

^-    =  ^^-..-^    -  ^^  <...^  ^3j 

Treating  the  booster,  M^,  as  a  free  body,  the  sum  of  the 
vertical  forces  acting  on  it  are 

TV  =   ^«>=ff^  ~P^,>f ^  -  ^ - ^  ^  -  z*^^*  -  /^^ 
and  therefore 

Furthernnore ,  the  sum  of  the  horizontal  forces  acting  on  M^  are 

and  therefore 

4.      --    ^^^^^y      ^^^-=f^     -    ^^p     -    ^/^  (5) 

Finally,  the  sum  of  the  moments  acting  on  M^  are 
or 

if^  ^^^fi  =  /-^  ''"Y^  -^A  ^r^f,  -  j^fA^A)  (6) 

Substituting  the  results  of  rlquetlon  (1)  into  Equation  (4)  to 
eliminate  f, 

Similarly,  substituting  the  results  of  Equation  (2)  into  Equation  (5) 
to  eliminate  f;^ 
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Substituting  the  results  of  lilquations  (4)  and  (5)  into  Equation  (3) 
to  eliminate  f ,    and  f^ 

(9) 

Similarly,  substituting  the  results  of  Equations  (4)  and  (5)  into 
Equation  (6)  to  eliminate  f,    and  f^ 

-^/-f'^;  (10) 

Assuming  the  angles  q^  and  q^  sufficiently  small  such  that 
terms  of  order  q*  may  be  neglected,  the  following  approximations 
can  be  made 

i/'r?^.  <b.   ^,  and  ^a^  ^.  •Ji  / 

Therefore,  Equations  (7),  (8),  (9).  and  (10)  may  be  rewritten  In 
sinnpler  form.     Thus 
From  Equation  (7) 

From  liquation  (8) 

From  Equation  (9) 

(13) 
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From  Equation  (10) 

The  following  approximate  relations  between  the  principal 
coordinates  q     ,  q     »  %  »  ^^  <ly  »  *nd  the  extraneous  coordinates 


and 


fj-     ^    /^    '^  /^^%i   '     '/^         "'/^Z 


(15) 


^^^.r^^-^'^^  <^^) 


^      -    /T    -^^^J^fs 


(17) 


Taking  the  first  and  second  time  derivatives  of  Equation  (15) 
Taking  the  first  and  second  time  derivatives  of  Equation  (16) 


raking  the  first  and  second  time  derivatives  of  ::^quation  (17) 
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Taking  the  first  and  second  time  derivatives  of  ilquation  (18) 
^     -     ^    ^   ,^-^>       rj-J^)  (25) 


and 


fs     '  r^    -^f^f^^-^^  (26) 

IDquation  (11)  can  now  be  rewritten  In  terms  of  the  principal 

coordinates  only  by  substituting  the  results  of  i£quatlons  (22),  (25) 

and  (26)  for  the  extraneous  coordinates 

Note  that  Elquatlon  (27)  Is  Identical  with  Equation  (27)  of  Part  II  of 
this  paper. 

In  a  similar  manner.  Equation  (12)  can  be  rewritten  In  terms 
of  the  principal  coordinates  by  substituting  the  results  of  Equations 
(20),  (23),  and  (24)  for  the  extraneous  coordinates 

which  Is  Identical  with  Equation  (23)  of  Part  II  of  this  paper. 

Substitution  of  Equations  (23),  (24),  (25),  and  (26)  In  Equa- 
tion (13)  yields 
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which  ia  Identical  with  H^quation  (44)  of  Part  II  of  this  paper. 

Finally,  aubstitution  of  Equations  (23),  (24).  (25),  and  (26) 
in  i£quation  (14)  yields 


which  is  identical  with  Equation  (59)  of  Part  II  of  this  paper. 
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TABLii:  OF  SYMBOLS 

-  -f 

^  -f 

D  -  Jet  damping  force 

F  -  Thrust  force 

g  -  Acceleration  of  gravity 

I  -  Moment  of  inertia 

I  -  Specific  Impulse 

k  -  Radius  of  gyration 

M  -  Mass 

M  -  Mass  rate  of  flow 

q  -  Lagrangian  generalized  space  coordinate 

Q  -  Dimensionless  Lagrangian  generalized  space  coordinate 

^  -  Lagrangian  generalized  force 

T  -  Kinetic  energy 
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